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ABSTRACT 

We present a generalized theory of microphase separation for charged-neutral diblock copoly- 
mer melt. Stability limit of the disordered phase for salt-free melt has been calculated us- 
ing Random Phase Approximation (RPA) and self-consistent field theory (SCFT). Explicit 
analytical free energy expressions for different classical ordered microstructures (lamellar, 
cylinder and sphere) are presented. We demonstrate that chemical mismatch required for 
the onset of microphase separation {x*N) in charged-neutral diblock melt is higher and the 
period of ordered microstructures is lower than those for the corresponding neutral-neutral 
diblock system. Theoretical predictions on the period of ordered structures in terms of 
Coulomb electrostatic interaction strength, chain length, block length, and the chemical 
mismatch between blocks are presented. SCFT has been used to go beyond the stability 
limit, where electrostatic potential and charge distribution are calculated self-consistently. 
Stability limits calculated using RPA are in perfect agreement with the corresponding SCFT 
calculations. Limiting laws for stability limit and the period of ordered structures are pre- 
sented and comparisons are made with an earlier theory. Also, transition boundaries between 
different morphologies have been investigated. 
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I. INTRODUCTION 



Science behind the complex behavior of amphiphihc systems continues to be of in- 
terest to scientific community. A great deal of theoretical^'^'^'^'^'^'^'^ and experimental 
effort8^>^>^>^'^>^ have been made to study this behavior. Especially, self-assembly of 
amphiphiles''^^ is of significant importance in understanding many biological systems. An 
amphiphilic diblock copolymer system has applications such as encapsulationp^i^ and drug 
deliveryi^, which are dependent on self-assembly of macromolecules. 

Self-assembly of neutral block copolymers in concentrated regimes, has already been 
studied extensively in the last three decades. Seminal work in developing theory for mi- 
crodomains in diblock copolymer melt was done by Helfand^^ et al , where unit cell ap- 
proximation was used to calculate the properties for sharp interfaces (strong segregation 
limit (SSL)). Later on, Leibler— calculated morphology diagrams in weak segregation limit 
(WSL) using the random phase approximation (RPA), when the system is on the verge of 
transformation from disordered to ordered phases and interfaces between the ordered do- 
mains are diffuse. Few years later, Semenov^^ and Ohta^^ et al calculated morphology 
diagrams in SSL taking into account a sharp interface between the domains. Semenov used 
ground state dominance and Ohta et al extended the RPA method developed by Leibler 
to SSL. Afterwards, Muthukumar— et al and Matsen^ et al bridged the gap between 
WSL and SSL theories by using density functional theory (DFT) and SCFT respectively. 
To go beyond mean field, fluctuations of order parameter were included by Fredrickson^ et 
al , Olvera de la cruz^° and Muthukumar— . 

Although we have a sound understanding of neutral copolymers, our understanding of 
charged copolymers is inadequate. A number of researchers have tried to explore charged 
diblock systems. Some of these efforts are invested in studying dilute solutions (micelle 
regime )^i2^ and others have been carried out in the concentrated regime^i^i^i2ii2^. The fun- 
damental question is how Coulomb interactions affect the relative stabilities of ordered 
morphologies. First work in this direction was carried out by Marko and Rabin^ who ex- 
plored charged copolymers in both the melts and solutions. In their study, they presented 
the effect of degree of ionization and salt on critical parameters for weak segregation in melt 
and studied micelle behaviour also. But they did not consider any ordered microstructures 
for melts. Recently, SCFT for polyelectrolytic systems has been developed by Shi^^ et al 



and Wang2^ et al . 

In this paper, we have considered microphase separation in diblock copolymer melts, 
formed out of polyelectrolyte [A) block and neutral {B) block. RPA is presented in Sec. Ill Al 
and SCFT equations are briefly presented in Sec. IIIBI Calculated results and conclusions 
are presented in Sec. [1111 and HVl respectively. 

II. THEORY 

We consider a system of n block copolymer chains, each containing a total of segments 
of two species {A and B) with / as the fraction of the block A. Number of segments in blocks 
A and B is represented by Na and Nb, respectively, so that (/ = Na/N,N = Na + Nb). 
Block A is taken to be polyelectrolytic (negatively charged) and we assume that there are Uc 
counterions released by the charged block. In addition to this, there are ions of species 
7 coming from the added salt (in total volume Q ) so that the whole system is globally 
electroneutral. Let Zi be the valency of the i*^ charged species. Subscripts A, B, c, + and 
— are used to represent monomer A, B, counterion from block A, positive and negative salt 
ions, respectively. The degree of ionization of the A block per chain is taken to be a so that 
each of the faN segments of A block chain carry a charge of cZa, where e is the electronic 
charge. So, there are A^(l — fa) uncharged segments in a chain. Recently, it has been shown 
that specific charge distributions along the backbone play a significant role in the physics of 
polyelectrolytes^. In our current study, we consider smeared charge distribution i.e. each 
segment of block A has charge cuZa- We represent copolymer chain as a continuous curve 
of length Nl, where / is the Kuhn segment length. For the treatment shown below, we 
assume that volume occupied by each A and B monomer is the same (= = 1/ poiPo being 
the bulk density) and that the system is incompressible so that Vt = nNl^ . We use an arc 
length variable Sj to represent any segment along the backbone of j*^ chain. In this notation, 
A block in j*^ chain is represented by Sj (0 < Sj < NaI)- Also, the position vector for a 
particular segment is represented by Rj{sj). For this system, Helmholtz free energy F can 
be written within Edward's formalism by: 
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Here, kBT is Boltzmann constant times the absolute temperature. 
^yi(^)) ^bb(^) and V^B(r) are the interaction energies between segments of different 
types separated by distance r =| r | so that 

VAA{r) = WAAS{r) + ^rrrfr- 2 

VBB(r) = WBB^ir) (3) 
Kii?(r) = WAB5(r) (4) 

where waa-, wbb and wab are the intersegment excluded volumes arising from the short range 
interactions, 6{r) is Dirac delta function and e is position independent effective dielectric 
constant of the medium ( in units of Aneo where eo is the permittivity of vacuum). In 
writing VaaI'"), it is assumed that total charge of any chain is uniformly distributed along 
A block (smeared charge distribution). Interactions between polymer segments and small 
charged molecules (counterions and colons) are accounted for by treating small ions as point 
charges so that they have zero excluded volume and are purely electrostatic in nature. These 
interactions, represented by VAm{r), VBm{r) and Knp(r) in Eq. ([!]), are given by 

VBmir) = (6) 
Vr^Ar) - '-^\ (7) 

The Dirac delta function involving pAi'f') and Pb{i^) enforce the incompressibility of the 
system at all locations, where PaI^) and Pb{t) are microscopic monomer number densities 



defined as 
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A. WSL - Random Phase Approximation 

1. Stability Limit 

Near the stability limit of the melt, densities of all components in the inhomogeneous 
phase deviate slightly from their average values in the homogeneous phase. So, the free 
energy of the inhomogenous phase can be obtained by expanding the corresponding free 
energy expression for the homogeneous phase about the average densities. Neglecting cubic 
and higher order terms in densities, the free energy of the inhomogeneous phase is obtained 
in terms of the order parameter 0(r) (see Appendix A) and has the form: 

F = Fo + 6F{(f)) + 0{<p'), (10) 

where the order parameter is the same as that used for neutral block copolymer melts, 
given by 

(P{r) =<6pA{r) > = < pAir) > - f po. (11) 

Fq is the free energy of the homogeneous phase and is given by a Flory-type equation: 

Fn ^ , , . 



k^T ^ n.ln(n.)- — + W(l-/) (12) 

i=c,+,— 

In writing Fq, linear and constant terms have been ignored. In Eq. ( [T2l) . xab and k,~^ = td 
are Flory's chi parameter and the Debye screening length, respectively, is given by 

and Ib = e^/ e/c^T stands for the Bjerrum length (in units of Aireo). In Eq. ( [T2ll . the first term 
corresponds to the translational entropy of small ions, term accounts for counterions and 
colons correlation effect and xab term gives the contribution of chemical mismatch between 
blocks. 



The second degree term in defines the structure factor {S{k)) and is given by: 

SFi<p) = -j j^,S-\k)m<f>i-k) (14) 

S-\k) = So\k) + S^\k) (15) 

S,\k) = Q{x)-2xabI^ (16) 

^ PoN {g{f, x)g{l - f,x) - [g{l,x) - g{f, x) - g{l - /, x)Y/A} ^ ^ 



- Qx + Kmp ^^^^ 



In Eq. ( [151), SQ^{k) is the contribution due to short range excluded volume interactions 
and Si^{k) is due to the long range Coulomb interactions present in the system. g{f,x) is 
the Debye function given by 



2(/x + e-f^ - 1) fcW 

g{f, x) = , X = — — = k R (19) 
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Note that Eqs. ( [14] - [18]) are the same as in Ref .- 



2. Limiting Laws 

To understand the qualitative behaviour of the system at the stability limit, a scaling 
analysis is presented here for two limiting cases. In order to carry out the scaling analysis, 
the complicated equation for SQ^{k) is approximated by an expression, used previously for 
neutral block copolymersli, which reproduces S^^^k) with about 5 % accuracy. Using the 
approximate expression, Eq. ( [H]) becomes: 

«(/) (21) 
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^(/) = Npp{i-fy ^^^^ 

C{Ib) = AnlsZla^poN (23) 

s(f) 

X = 2xabNpoI' - (24) 

Here, s{f) is a parameter used to reproduce SQ^{k). At the stability limit, S{k) diverges 
at the wave-vector k = k* and the second degree term 5F in free energy vanishes aX k = 
k* axidxAsN = Xab-^- ^^^e of neutral block copolymers (i.e. C{Ib) = 0), the 



divergence corresponds to k* — [A/By/^ ~ N ^1"^ so that the period D — li^jk'^ ~ N'^l'^. 



Hence, vanishes at %* = 2\J AB so that x\b^ is given by expression 

\XaB^^ ) neutral 



(25) 
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Now, for charged block copolymers, consider the two limiting cases: (i) k"^ <^ and (ii) 
k'^ ^ K^. In terms of the Debye screening length , these cases correspond to -C -D and 
ro ^ D, respectively. 

Case (i): td ^ D 
In this regime, the structure factor becomes 

S-\k)c.Bk-^^^-[x-^)) (26) 

The maxima of the structure factor corresponds to A;* = [A/BY/^ ~ N~'^/'^ so that the 
period D — 21: /k* ~ N^^"^ and is independent of degree of ionization of charged block. 
So, charged block copolymer behaves like neutral copolymer. This is because the electrostatic 
interactions are short ranged in this regime. But there is a remarkable effect of small ions on 
Xab-^ and in fact, it is found that charged block copolymers have to have higher Xab^ ^ 
compared to its neutral analog for undergoing the microphase separation. For salty systems, 

{X*AB^) charged ^ bC*AB^) neutral ^2 (^''') 

For salt-free incompressible system where — AnlBfaZ^/l^, this expression simplifies to 

1 fZAVaN 



(XaB^) charged " (X AB ^) neutral + 2 ) ~~f ^^^^ 

Physically, this means that homogeneous phase in charged copolymer enjoys larger param- 
eter space as compared to its neutral analog. 

Case (ii): rn ^ D 

In this regime, the structure factor becomes 

S-\k) = Bk' + ^^-X (29) 



The maxima of the structure factor corresponds to 



(30) 



so that X* = 2>/ (A + C)B. Plugging in expressions for A, B and C 
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(32) 



From the expression for it is clear that the period decreases with an increase in the 
degree of ionization in this limiting case. Comparing Xab-^ charged and neutral block 
copolymer cases (Eq. ( [32]) and Eq. ( [25!) ). it can be inferred that in this limiting regime 
also, {x\bN) charged > (x^bN) neutral ^nd iucrcascs with an increase in a. 

Further, it can be shown that the correlation effect of ions (Debye-Hiickel theory) is weak 
as long as kZb ^ 1. This means that all of the above limiting laws for salt-free systems are 
valid as long as kIb <^ 1. 



3. Numerical Calculations for Stability Limit 

From here onwards, we consider the salt-free melt. To calculate the stability limit for 
salt-free melts, inverse temperature dependence of xab and Ib is clubbed together by the 
introduction of a parameter (called reduced temperature^) defined as 



Having written structure factors in terms of t, the stability limit is calculated using 

U=.^.=** = (34) 

OK 

S-\k) \k=k*,t=t* = (35) 



where 5* ^{k) is given in Eqs. ([I5]- [TSl) . Solving these equations for t* 



where 



, -L + y/L^ - APR , , 

= -1^ > 0, (36) 



P = 607rxM(x*) (37) 

L = lOn fNZ^A{x*) - 6x* + lOnZ^aN (38) 

R = -fNZl (39) 

A{x*) = Q{x*) (40) 



In writing these equations, we have taken / = 1. Function Q appearing in these equations 
has aheady been defined in Eq. ([I7|). Using Eq. ([Ml), the wavevector at the stabihty hmit 
is given by solving the equation 

SQjx) GZlaNt* 

First, Eq. ( I^Tl) is solved for x* and then t* is calculated using Eq. ( 1361) . Effect of a and 
on the stability limit is shown in Figs. [H - O obtained by using Eq. ( HTl) for different 
values of N is shown in Fig. [3l 

4- Ordered Structures 

To derive free energy expressions for different ordered structures, we employ the method 
used by Leibler— . Following Leibler, we expand the free energy expression in terms of the 
order parameter up to fourth order. Taking advantage of the fact that in WSL (near stability 
limit), important fluctuations in polymer densities are those with the wavevector k = k*, 
we approximate the order parameter by a sum of plane waves, each having the wavevector 
k = k*. Using this expression for the order parameter, free energy density becomes 

SFn = = 2Nl'{Xs - X)<Pl - Cn€ + VnCt>'n (42) 

2Xsl' = g(x^) + g^,/^ and x = Xab (43) 

where the value of n corresponds to the morphology being studied. To be speciflc, n = 1,3 
and 6 correspond to lamellar, hexagonally close packed (HCP) cylinder and body centred 
cubic (BCC) spherical morphology, respectively. Functions Fs, r4 and coefficients Cn, Vn were 
calculated by Leibler— ( summarized in Table [T]). Further, we have adopted the notation 
used in Ref.— for the arguments of the function The coefficients have the property that 
Vi < Vs < V6 for all / and specifically, for / = 1/2, ^„ = 0. For all the calculations presented 
in this paper, these coefficients are evaluated at k = k*. 

5. Transition Boundaries 

By following the Leibler's procedure^, the disorder-order transition (DOT) and the order- 
order transitions (OOT) are studied. Minimizing the free energy density (Eq. ( W2\i ) with 



respect to the order parameter, equilibrium order parameter and free energy densities are 
obtained. Results of these minimizations are presented in Table [III where 7„ is given by 

.1/2 



In 



1 6^^" ( \ AT 



(44) 



In order to determine the morphology that evolves at DOT, the free energy density is 
equated to zero so that XnN at DOT is found to be 

Xr.N = XsN-f- (45) 
8r]n 

Using the coefficients (n and rjn, it can be shown that BCC {n = 6) gives the lowest value 
for XnN (or highest value of t). So, the morphology that appears first is BCC. Writing Eq. 
( l45l) for DOT in terms of the reduced temperature t, the DOT boundary is given by Eq. 
( En]) where A{x*) is now given by 

Subscript 6 implies that the morphology is sphere {n = 6). Similarly, order-order transition 
boundaries are calculated by equating free energy densities for different morphologies. In 
general, all the transition boundaries (stability limit, DOT and order-order transitions) are 
calculated using Eq. ( 1361) . where only A{x*) varies. Mathematical conditions and values of 
A{x*) for different transition boundaries are summarized in Table IIIII Function y appearing 
in Table IIIII is the solution of Eq. (V-35) in Leibler's work^^. Solving these sets of equations, 
the morphology diagram can be constructed as discussed in Sec. IIIII 



B. Self- Consistent Field Theory (SCFT) 

Although RPA gives a valuable insight into the physics of the problem, it is only a linear 
response treatment. Strictly, this treatment is valid close to the stability limit of homoge- 
neous phase but far from the limit, RPA calculations are not quantitatively correct^i^^i^. To 
go far away from the stability limit, SCFT has been used extensively in the literature^'^'^. 
Using standard methods^^, self-consistent equations are obtained under the saddle point 
approximation (see Appendix B). 

This numerical technique leads to coupling of full non-linear Poisson-Boltzmann equation 
with standard modified diffusion equation for the polymer chains. We have solved these 



sets of equations using an efficient spectral technique^!. While solving these equations, 
experimentally found inverse dependence of Flory's x parameter has been exploited by using 
reduced temperature t (Eq. ([33])). We have studied the effect of a on Xab^ compared 
with the corresponding RPA calculations. Also, the effect of degree of segregation on the 
period of lamellar morphology is studied (Fig. Hj). Monomer densities, counterion densities 
and electrostatic potential obtained from SCFT calculations are shown in Figs. O - El 
respectively.. 

III. RESULTS 

In the previous section , we have provided the necessary equations to describe the mi- 
crophase separation. Here, we present results for salt-free charged-neutral diblock melts, by 
solving the above equations. 

A. Stability limit - RPA results 

In Figs. [Hand [21 we have drawn the stability limits for charge/neutral block copolymer 
salt-free melt at different degrees of ionizations for N = 1000 and = 10, 000, respectively. 
The critical value of the reduced temperature required to induce microphase separation de- 
creases with an increase in the degree of ionization. This is in qualitative agreement with 
limiting laws presented in Sec. Ill Al and the already established concept that the effective 
Flory's x parameter decreases with an increase in degree of ionization for polyelectrolytes^^. 
Unlike the neutral copolymers, x ^^^^ ^ for polyelectrolytic diblock copolymers are inde- 
pendent parameters that govern the phase behaviour. Also, the stability limit depends on 
fraction of charged block (/) in an unsymmetric fashion. It is to be noted that these results 
are in agreement with the results reported by Marko and RabinI'. In Ref.-, temperature de- 
pendence of X parameter was not taken into account and critical parameters were calculated 
by choosing a fixed value oils /I- The method of calculations used by Marko and Rabin was 
similar to the one presented in Appendix A. 



B. Period of lamellar phase (/ = 1/2) 



In RPA, near the stability limit, the period of an ordered structure is approximated by 
D = 27r/A; ~ 27r/A;*. It is well known that the mean field theory^^ for neutral block copolymer 
predicts Xab-^ — 10.495 and x* = 3.7852 at / = 1/2. Also, the period shows 1/2 power 
law dependence on in WSL (i.e. x = constant) as long as the wavevector dependence 
of higher order terms in free energy expression is suppressed^i^i^. Physically, this means 
that block copolymer chains are obeying the Gaussian statistics for chain conformations. 
Experimentally, there are deviations from this power law because of chain stretching^^. As 
shown in Fig. [3l same power law dependence is obtained using RPA for polyelectrolytic 
block copolymers when A^ is large, but the period {D ~ ^N/x*) for a given A^ is smaller 
than that for an equivalent neutral copolymer system. This effect has been seen by other 
researchers also^. The decrease in period with increase in degree of ionization, is explained 
in Ref.- by an argument that counterions need to be rearranged on microphase segregation 
and entropy loss is lower if the length scale of fluctuations for these counterions is smaller. 
It is to be stressed that lowering of D with a is not purely entropic effect. This effect is an 
outcome of electrostatic screening due to counterions and hence, includes both energetic as 
well as entropic contributions. 

At present, we are not aware of any experimental data on the period of charged- neutral 
block copolymer. Nevertheless we expect the polyelectrolyte chains to be non-Gaussian and 
D to deviate strongly from N^^"^ power law in the case of charged systems. SCFT has 
been quite successful in predicting the period of ordered structures for neutral copolymers. 
Expecting that SCFT results are valid for weakly charged polyelectrolyte copolymers, results 
obtained from SCFT calculations are plotted in Fig. IHfor lamellar phase (/ = 1/2, N = 
1000). Lowermost point in the plots of Fig. H] corresponds to Xab-^- comparing Fig. 
Hlwith Figs. [Hand El it is clear that the RPA calculations for Xab-^ good 
agreement with the corresponding SCFT calculations. Analogous to neutral copolymers, 
it is found that N^^'^ power law is not valid for ordered microstructures. In addition the 
quahtative feature that the period of ordered structures is lower than its neutral analog, is 
clearly seen in these plots. 



C. Counterion distribution 



The RPA calculations do not provide counterion distributions. On the other hand, SCFT 
allows calculations of counterion densities and potential self-consistently. Figs. [5] - [7] 
show monomer densities for charged block (A), counterion densities, and the electrostatic 
potential, respectively. The onset of microphase separation leads to creation of monomer 
and counterion density waves (Figs. and E]) because of the incompatibility between the 
blocks and coupling between charged monomer (A) and counterions, respectively. From 
these plots, it can be inferred that in the strong segregation limit {xN oo), all the 
counterions are confined to the charged domains. One of the effects of these density waves 
in the lamellar phases, is the presence of a potential difference between the charged and 
neutral domains (Fig. [7]) whose magnitude increases with the degree of segregation. For 
the weakly charged diblock copolymer system studied here, total local charge density is 
close to zero and it is hard to determine the shape of the charge density wave for the system 
due to the possible numerical errors. To verify the observation that the effective degree 
of segregation is reduced because of the electrostatic interactions (RPA calculations), we 
have plotted monomer densities for neutral and charged block copolymer at the same 
(Fig. [8]). These plots clearly confirm that effective degree of segregation is lower for charged 
copolymer melt and counterions have a tendency to drive the system towards homogeneous 
phase. 



D. Morphology diagram for charged- neutral diblock copolymer 

The calculation of morphology diagrams using SCFT is a computationally intensive task 
because of the vast parameter space for polyelectrolytic systems. To get an idea about 
transition boundaries, we have used RPA calculations presented in section III A| assuming 
that only the classical morphologies^ compete in charged-neutral diblock copolymer sys- 
tems. Figs. M and [TU] show the calculated morphology diagrams for different a and A^. We 
observe that DOT and OOT boundaries are strongly dependent on a and A^. The temper- 
ature of occurrence of DOT decreases with an increase in a and increases with an increase 
in A^ (analogous to the shift of stability limit - Figs. [H [2]) . Furthermore, these transition 
boundaries for DOT and OOT are highly asymmetric with respect to /. 



IV. CONCLUSIONS 



We have addressed the microphase separation in charged-neutral diblock copolymer melts 
in the weak segregation limit, by using the RPA and SCFT methods. We have shown in 
Sections III Al and IIIII that the critical value of the x parameter for microphase separation 
is higher for charged copolymers and of concentration modulation is smaller in comparison 
with neutral copolymers. From morphology diagrams, it can be seen that the parameter 
space for ordered microstructures is reduced when degree of ionization of charged block is 
increased. In other words, charging a block stabilizes the homogeneous phase. 

The SCFT results show that the counterions partition themselves preferentially within 
the charged domains. This leads to creation of a potential difference between charged and 
neutral domains. This process of partitioning is unfavorable both entropically and energet- 
ically. Hence, the length scale of these partitioning is lower when there are more number of 
counterions to be partitioned for the same number of monomers. 

Finally, we summarize the assumptions in obtaining the above results. We have taken the 
counterions to be point charges. Our treatment can be readily extended to counterions with 
finite size by modifying the incompressibility condition and incorporating excluded volume 
interactions. The Kuhn segment lengths for the neutral and charged blocks are taken to be 
same. It has been shown that conformational asymmetry has an effect on the order-order 
transition boundaries for neutral copolymers'^. Analogously, there will be an effect on our 
system as well. In the case of polyelectrolytes, electrostatic interactions cause stiffening of 
the chain so that the effective segment length^ depends on various factors such as k, a etc. 
in a complicated manner. By assuming that the segment length for charged and neutral 
blocks to be the same, we implicitly assume that the charged block copolymer under consid- 
eration is conformationally symmetric and is weakly charged, so that the difference between 
the effective and bare segment lengths is negligible . Another important assumption in the 
present theory is that the position/concentration dependence of the dielectric constant e is 
suppressed . Recently, the effect of dielectric constants of individual components in a multi- 
component polyelectrolyte^ solution has been presented. In this work, dielectric constant of 
a component was taken to be linearly dependent on the concentration of the component. In 
principle, dielectric constant depends on the concentration of ions in a complex manner—. 
At present, there is no satisfactory well-established model for the dependence of microscopic 



dielectric constant on macroscopic density. So, we model dielectric constant (e) appearing in 
the expression for the Bjerrum length as the effective dielectric constant for the mixture of 
A, B monomers and counterions. Position dependence of dielectric constant will definitely 
play an important role in the strong segregation limit. However, for melts in WSL, the 
average value of dielectric constant can be taken as a constant. Further, we have consid- 
ered only the lamellar, cylindrical and spherical morphologies as the competing structures. 
Extensions of the present theory for other morphologies are in progress. 

At present we are unaware of any systematic experimental study on charged-neutral 
copolymer. We hope that the present theoretical work will instigate experimental work on 
charged block copolymers in the concentrated regime. 
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APPENDIX A : Integration over positions of small ions 



Here, we present how to compute the integral over positions of small ions. If microsopic 
charge density at any point in space is defined as 
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where Pyi(r) is given by Equation [8] and charge density of small ions is defined as 

Pm{r) = J2z,^6{r-r,). (A-2) 
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Using these definitions of densities. Equation [T] can be written in the form 
exp 
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Here, xab is Flory's chi parameter, which accounts for chemical mismatch between blocks 
and is defined as 

XabI =wab ^ (A-4) 

Terms corresponding to excluded volume interactions can be treated in the usual way^. 
Here we present the treatment of electrostatic terms. We introduce field variables to go from 



particles to fields by defining 0p(r) 



PO 



where p = A, B, c, +, — (i.e. p corresponds to 



all components in system i.e. y4-block,i?-block,counterions and colons). Before proceeding, 
we introduce the identity 
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Here, i is purely imaginary number and to avoid confusion, we are going to use the notation 
iwp{r) Wp{r) for all purely imaginary fields. In this notation. Equation IA-31 can be 
written in the form 



exp 



l[D[(Pp{r)]D[wp{r)]D[r]{r)]exp 



H 



(A-7) 



exp 



H 



T 



7 "7 



exp 



^Vo / drxAB(t)A{r)(t)B{r) 



Pi 



+Po j dr\S^Wp{r)(l)p{r) +ri{r) (^^(r) + ^^(r) - 1) 



dr I df 



/ pe{r)pe{r') 
eksTlr — r'\ 

QIbQ^'QI^Q-- 

(A-8) 



where ir]{r) r]{r) is the field introduced to enforce the incompressibility condition at all 
points and 



Pe[r) = e 



ZAfa(j)A{r) + ^ 



m=c,+, 



(A-9) 



Qab = - I D[R]exp 



3 , fdR(s] 

'21 L ''[-^ 



1 



1 



ds WAiRi-s)) — - / dswBiRi-s)) 

' JnaI 

(A-10) 



Nl 



Here, subscript i{= c, +, — ) corresponds to all kinds of small ions in the system. Q is the 
total volume and r is a normalization constant given by 



D[R]exp 



Nl 

I ds 

21 



dRjs) 
ds 



(A-12) 



Integrals over Wp can't be calculated exactly. So, in order to proceed further, integrals over 
Wi{r) are approximated by the maximum value of the integrand (saddle point approxima- 
tion). Maximization of the integrand with respect to Wi{r) gives 



rii 



Plugging Wi back into Eq. ( I A- 71 
F 



exp [— u'i(r)] ^ Wi{r) 



-In 



Hi 



In [0i(r)] 



exp 



exp 



\{D[(t)p{r)]D[wA{r)]D[wB{r)]D[r]{r)]ex^ 
p 



(A-13) 



(A-14) 



exp 



7 '"J 

I dr I dr 



-I po / drxAB(l)Air)(l)Bir) 



, pe{r)pe{r') 



eksTlr — r' 



- + Po / drwA{r)(t)A{r) 



+Po / drwB{r)(f)B{r) - po dr ^ 0j(r) In [0i(r)] 

^ i=c,+ 



Qab- 



(A-15) 



+Po J drr]{r){(j)Air) + (j)Bir) - 1} 

In order to calculate the integral over 0j, we consider fluctuations of charge densities about 
the disordered phase- . For weak fluctuations, the entropic term 0jln0j can be expanded in 
powers of these fluctuations. Keeping the leading terms (neglecting cubic and higher order 
terms), integrals over 0j become Gaussian and can be easily calculated. So, by writing 

0i(r) =< (Pi > +6(Pi{r) (A-16) 

, and using J dr6(f)i{r) = in combination with the global electroneutrality in the homoge- 
neous phase, integrals over (pi (i.e. small ions) become 



J exp 



m=l 
3 



Pi 



d^k 



(27r)3 

+2 vt''^aZAZm5(pA{k)5(Pm{k) 



3 3 
, m=l q=l 



5, 



mq 



+ vt''^Z^ZA5(Pq{k) 



PO < (pm> 



m=l 



exp 



Pi 



d^k 
(2^ 



V^'^a^Zl6(P\{k) 



(A-17) 



where V, 



(0) 



— -p^, 6mp is the Kronecker delta and the indices 1,2,3 correspond to counterion 
from polymer, positive (+) and negative (-) salt ions, respectively. By calculating the 



Gaussian integrals, the value of / is found to be 



(2vr) 



P0VC1C2C3 



exp 



exp 



'2 



In 



Po ^2 72 
-—a Zj A 



(0) 



(2^)' 1 + T/fco 



where Cr, 



and Co is given by 



Co 



- y. 



7 ^ T) 



m=l,2,3 



B 



Now, using the relation 



X 



In 



1 + 



a 



a 



x^ 



TT 



-a 



we get 



2" 



f2 [ d^k AttIbCq 
(27r)3 P 



(A-18) 



(A-19) 



(A-20) 



(A-21) 



127r ^ ^ 2 

Here, the last term is divergent. But the divergence is for k being large (ultraviolet di- 
vergence) and we are interested on a length scale corresponding to small k. So, this term 
is neglected. Note that the same result can be obtained by employing the central limit 
theoren>2i. Now, electrostatic terms are decoupled from the rest of the terms and remaining 
treatment is the same as that done for the corresponding neutral copolymer—. 



APPENDIX B : Self- Consistent Field Theory 

Instead of expanding about average densities, the introduction of another field corre- 
sponding to charge density^^^ in Eq. ( IA-71) . leads to the following self consistent equations 
for the salt-free melt after taking / = 1 



XABNcpBir) = WA{r)+rj{r) 
XABN(f)A{r) = WB{r)+r]{r) 
4>A{r) + (pB^r) = 1 



(B-1) 
(B-2) 
(B-3) 



(t>A{r) 
(t>B{r) 

VrV'(r) 

dq{r, s) 
ds 

dq*{r, t) 
di 



ZAfa 

ZcQc 



exp [-Zcip{r)] 



Qj^ dsq{r, s)q*{r,l-s) 

Jdrq(r,l) 
Q ds q{r, s)q*{r, 1 — s) 

Jdrq(r, 1) 
-AttIb [Zc(f)c{r) + ZAa(f)A{r)] 



[f V? - {ZAaNij{r) + w^(r)}] q{r, s) s < f 

[^^l-WB{r)]q{r,s) s>f 

[f V? - {ZAaNiPir) + w^(r)}] g*(r, t) t > (1 - /) 



(B-4) 
(B-5) 

(B-6) 
(B-7) 

(B-8) 



(B-9) 



^<(1-/) 



These equations are to be solved with initial conditions q{r, 0) = 1, q*{r, 0) = 1 and the free 
energy expression for the salt-free melt (per chain) becomes 

F If, N 



nksT 



1 

- iln 



drxABN(f)A{r)(f)B{r) 



Q 



AB 



n 



ZAfaN 



In 



Qc 



j dr\\/ri}{r)\'^ + J drri{r 



rir 



ZAfaN 



(InQ + 1) + (Inr + 1) } (B-10) 



Here, we have used the notation ieip{r) — > ip{r) for purely imaginary electrostatic potential 
and iwp{r) — > Wp{r) for all the purely imaginary fields. Expressions for Qc, Qab and r are 
given in Appendix A. The free energy of the homogeneous phase (where all the densities 
and fields are constant) is given by 



nkBT 



XABNfil - /) + -^faN 



In 



Za. ' 

— -^fOL 

Z/ 



- 1 



(B-11) 



REFERENCES 



1 A.R. Khokhlov, J.Phys.A: Math. Gen.,13, 979 (1980). 

^ V.Yu. Borue and I.Ya. Erukhimovich, Macromolecules, 21, 3240 (1988). 
3 (a) J.F. Marko and Y. Rabin , Macromolecules, 24, 2134 (1991); (b) ibid. 25, 1503 (1992); 
^ M. Benmouna, T.A. Vilgis and J. Frangois, Makromol. Chem., Theory Simul., 1, 3 (1992). 
^ M. Benmouna and Y. Bouayed, Macromolecules, 25, 5318 (1992). 

^ E.E. Dormidontova, I. Ya. Erukhimovich and A.R. Khokhlov, Macromol. Theory and Simul., 
3, 661 (1994). 

G.A. Carri and M. Muthukumar, J. Chem. Phys., Ill, 1765 (1999). 

8 K. Ghosh, G.A. Carri and M. Muthukumar, J. Chem. Phys., 116, 5299 (2002). 

9 K. Nishida, K. Kaji and T. Kanaya, J. Chem. Phys., 114, 8671 (2001). 
K. Nishida, K. Kaji and T. Kanaya, J. Chem. Phys., 115, 8217 (2001). 
^ K. Yu and A. Eisenberg, Macromolecules, 31, 3509 (1998). 

2 S.A. Jenekhe and X.L. Chen, Science, 279, 1903 (1998). 

3 X.L. Chen and S.A. Jenekhe, Langmuir, 15, 8007 (1999). 

^ T. Nakanishi, S. Fukushima, K. Okamoto, M. Suzuki, Y. Matsumara, M. Yokoyama, T. Okano, 

Y. Sakurai and K. Kataoka, J. Control. Release, 74, 295 (2001). 
^ (a) E. Helfand, J. Chem. Phys. 62, 999 (1975);(b) Macromolecules ,8, 552 (1975);(c) ibid. 9, 879 

(1976);(d) ibid. 11, 960 (1978);(e) ibid. 13, 994 (1980). 
^ L. Leibler, Macromolecules, 13, 1602 (1980). 
^ There are few misprints in Ref.^^. Corrections are: 

(1) Equation (V-33) should read 

= [San + 1 + 2(3a„ + l)^/^]/3{an - 1) 

(2) Equation (C-6) should read 

54(/, h) = {1- exp[-(l - /)x]}{l + eM-hfx)/{h - 1) - /i exp(-/x)/(/i - l)}/hx^ 

(3) In Appendix C, functions gi{f,h) and fi{f,h) where i = 2,3,4, are defined and these 
functions have special cases for /i = 0, 1, 3, 4. In fact, these functions differ from general expres- 



sions(like Eq. (C-4)) only for ^ = 0, 1 and for all other values of h, general expressions are true. 

So, equations like (C-4a) and (C-4b) are valid only for = and h = l respectively. 
^8 (a) A.N. Semenov, Zh. Eksp. Teor. Fiz.,88, 1242 (1985); (b) Macromolecules, 22, 2849 (1989). 
^9 (a) T. Ohta and K. Kawasaki , Macromolecules, 19, 2621 (1986); (b) ihid. 21, 2972 (1988); (c) 

ihid. 23, 2413 (1990). 

2° (a) E.B. Zhulina and O.V. Borisov , Macromolecules, 35, 9191(2002); (b) O.V. Borisov and E.B. 

Zhulina , ihid., 35, 4472(2002); (c) ihid., 36, 10029(2003). 

21 N. Dan and M. Tirrell, Macromolecules, 26, 4310 (1993). 

22 N.P. Shusharina, LA. Nyrkova and A.R. Khokhlov, Macromolecules, 29, 3167 (1996). 

23 LA. Nyrkova and A.N. Semenov, Faraday Discuss. , 128, 113(2005). 

24 A.V. Kyrylyuk and J.G.E.M. Praaije, J. Chem. Phys., 121, 2806(2004). 

25 A.V. Kyrylyuk and J.G.E.M. Fraaijc, J. Clicm. Phys., 121, 9166(2004). 

26 M. Muthukumar, J. Chem. Phys., 86, 7230 (1987). 

27 M. Muthukumar , J. Chem. Phys., 105, 5183 (1996). 

28 M. Muthukumar, Macromolecules, 35, 9142 (2002). 

29 G.H. Fredrickson and E. Helfand, J. Chem. Phys.,87, 697 (1987). 
3° M. Olvera de la Cruz, Phys. Rev. Lett. ,67, 85 (1991). 

31 J. Melenkevitz and M. Muthukumar, Macromolecules,24, 4199 (1991). 

32 R.L. Lescanec and M. Muthukumar, Macromolecules ,26, 3908 (1993). 

33 M. Muthukumar, Macromolecules ,26, 5259 (1993). 

34 K. Almdal, J. H. Rosedale, F. Bates, G. D. Wignall and G. H. Fredrickson, Phys. Rev. Lett. ,65, 
1112 (1990). 

35 A. Shi and J Noolandi, Macromol. Theory Simul.,8, 214 (1999). 

36 (a) Q. Wang, T. Taniguchi and G.H. Fredrickson, J. Phys. Chem. B, 108, 6733 (2004); (b) ihid. 
109, 9855 (2005). 

37 (a) M.W. Matsen and M. Schick, Phys. Rev. Lett.72, 2660 (1994); (b) M.W. Matsen and F.S. 
Bates, Macromolecules,29, 1091 (1996); (c) M.W. Matsen, J. Phys.: Condens. Matter 14, R21 
(2002). 

38 (a) J. D. Vavasour and M. D. Whitmore, Macromolecules 25, 5477 (1992); (b) M.W. Matsen 
and F.S. Bates, J. Polym. Sci. B: Polym. Phys. 35, 945 (1997); 

39 G.H. Fredrickson, The Equilihrium Theory of Inhomogeneous Polymers (Oxford University, New 



York, 2006). 

C.J.F. Bottcher, Theory of Electric Polarization (Elsevier, Amsterdam, 1973). 



FIGURE CAPTION 



Fig. 1.: Effect of degree of ionization (a) on the stability limit in polyelectrolytic diblock 
melt: plots correspond to iV = 1000 and a = 0, 0.01, 0.02, 0.1. 

Fig. 2.: Effect of degree of ionization (a) on the stability limit in polyelectrolytic diblock 
melt: plots correspond to iV = 10, 000 and a = 0, 0.01, 0.02, 0.1. 

Fig. 3.: RPA Calculations - Effect of degree of polymerization (A^) and degree of ionization 
(a) on critical parameter x* for (/ = i) in WSL : a = 0, 0.01, 0.02, 0.1. 

Fig. 4.: SCFT Calculations - Effect of degree of segregation on period of lamellae (/ = 
1/2, N = 1000). Semenov's Strong Segregation Theory (SSST)^'^ which predicts 
D/ i^N^/'^l) = 2 (SxN/Stt'^)^^^ is also drawn for comparison purposes. 

Fig. 5.: Polyelectrolytic block copolymer lamellae (/ = 1/2, a = 0.01, = 1000) - 
monomer densities. 

Fig. 6.: Counterion distribution in lamellar phase (/ = 1/2, a = 0.01, A = 1000). 

Fig. 7.: Electrostatic potential in lamellar phase (/ = 1/2, a = 0.01, N = 1000). 

Fig. 8.: Reduction of effective chemical mismatch (/ = 1/2, N = 1000) - comparison be- 
tween monomer densities. 

Fig. 9.: RPA Calculations - Morphology diagram for polyelectrolytic diblock copolymer: 
= 1000 and a = for the topmost four boundaries, a = 0.01 for the middle four, 
a = 0.02 for the next set and a = 0.1 for the lowermost four boundaries. 

Fig. 10.: RPA Calculations - Morphology diagram for polyelectrolytic diblock copolymer: 
= 10, 000 and a = for the topmost four boundaries, a = 0.01 for the middle four, 
a = 0.02 for the next set and a = 0.1 for the lowermost four boundaries. 



Morphology 


Cn 




Lamellar 


Ci = 


r/i = fr4(o,o) 


Cylinder 


Cs = -(2/3V3)iVr3 


773 = i(r4(o,o) + 4r4(o,i)) 


Sphere 


Ce = -(4/3\/6)7Vr3 

+2r4(o,2) +4r4(i,2)) 


r?6 = i(r4(o,o) + 8r4(o,i) 



TABLE L Coefficients and r/„ calculated by Leibler— 



Morphology 


Order Parameter ((/>„) 


Free Energy Density {6Fn) 


Disorder - Order Transition 


Cn/ (2%) 





Lamellar 


Vix - Xs)N/vi 




Sphere, Cylinder 


Kn{l+ln)/{8Vn) 


27C^(l + 7n)'(l-37n)/(4096r?3) 



TABLE IL Equilibrium order parameters and free energy densities 



Transition Boundary 


Mathematical Conditions 


A{x*) 


Stability Limit 


Sk \k=k*,t=t*- u 
S~^{k) \k=k*,t=t* = 


Q{x*) 


Disorder - Order 


d{SF„) 1 

d(f>„ '<l>n=4>n 

5F„(<^"„,t„) = 


Q{x*) - ClI (4iV%) 


Sphere - Cylinder 


5Fq = SFs 


Q{x*) + 2y/N 


Cy]iud(\r - Laiiw^lLu; 


SF,, = (SFi 


9(2(^2 _i)/(32A'/;3) 



TABLE IIL Description of different transition boundaries 
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